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1.  Introduction 


Lower  leg  injury  is  a  major  injury  mode  indueed  by  underbody  blast  loading  assoeiated  with 
improvised  explosive  deviee  (lED)  attaeks  against  ground  fighting  vehieles.^"^  Blast -mitigating 
fioor  mats  are  regarded  as  an  effeetive  eountermeasure  for  injury  prevention,^’®  and  the 
methodology  for  their  seleetion  and  optimization  is  in  immediate  need.  However,  the  loading 
meehanies  in  the  lower  leg  with  the  existenee  of  the  fioor  mat  are  not  fully  understood.  Thus,  we 
deeided  to  elarify  the  meehanies  through  the  Finite  Element  Analysis  (FEA).  In  this  approaeh, 
we  seleeted  the  eombination  of  the  Hybrid  III  50*  Pereentile  anthropomorphie  test  deviee’ s 
(ATD’s)  lower  leg  and  the  butyl -rubber  fioor  mat  with  their  well-known  eharaeteristies. 

There  are  several  finite  element  (FE)  material  models  for  the  butyl  rubber,  and  it  is  neeessary  to 
determine  whieh  model  is  to  be  used  in  the  FEA.  In  this  report,  we  present  the  aeeuraey  of  FE- 
material  models  for  the  butyl  rubber  through  the  eomparison  between  the  FEA  and  eomponent 
loading  test. 


2.  Experiment  and  Analysis  Methods 


2,1  Butyl  Rubber 

The  butyl  rubber  used  in  this  study  was  manufaetured  by  Humaneties  Innovative  Solutions.  This 
is  the  same  material  used  in  the  neek  eomponents  of  the  Hybrid  III  50*  Pereentile  ATD.  The 
dimensions  of  a  bulk  sheet  were  8  inehes  x  8  inehes  x  20  millimeters  (mm),  and  test  speeimens 
were  eut  off  from  the  sheet  in  the  shape  of  eylinders  ((t)3  inehes  x  20  mm)  as  shown  in  Fig.  1 . 
Before  the  eomponent  loading  test  was  eondueted,  eaeh  speeimen  reeeived  eompression  loading 
up  to  40%  deformation  against  its  thiekness.  This  preloading  proeess  was  intended  to  initialize 
the  material  status  eonsidering  the  Mullins  effeet,  and  neither  permanent  deformation  nor 
fraeture  oeeurred  on  the  speeimens  during  this  proeess. 
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Fig.  1  Test  specimen  of  the  butyl  rubber 

2,2  Component  Loading  Test 

The  component  loading  test  was  conducted  on  the  drop  tester  with  the  finger-crusher  test  rig  in 
Fig.  2.  The  dimensions  of  the  test  rig  are  shown  in  Fig.  3.  In  this  test  setup,  the  lower  plate,  load 
cell,  and  base  were  rigidly  assembled  and  fixed  on  the  loading  table  of  the  drop  tester.  Thus,  the 
specimen  was  compressed  between  the  plates  in  the  drop  impact,  and  the  material’s  behavior  was 
known  from  the  applied  load  and  displacement  between  the  plates. 

The  test  was  conducted  on  the  drop  testers  (Lansmont,  P65  and  P45)  at  the  Adelphi  Laboratory 
Center,  U.S.  Army  Research  Laboratory  (ARL),  at  room  temperature.  The  loading  conditions  in 
the  test  are  listed  in  Table  1.  Test  specimens  were  reused  several  times  after  the  confirmation  of 
the  material’s  condition,  and  the  interval  between  each  test  was  kept  to  more  than  an  hour  for  the 
material’s  relaxation. 


Fig.  2  Front  view  of  tbe  finger-crusber  test  rig  on  tbe  loading  table  of  tbe 
drop  tester  equipped  with  a  specimen  and  sensors 
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381 


Mass  (Upper  Plate):  8.97kg 
Mass  (Lower  Plate):  4.48kg 

Bolted: 

Lower  Plate  /  Load  Cell  Interface 
Load  Cell  /  Base  Plate  Interface 


Top  View 


Upper 

Plate 


Lower  .  I-I  I  ^  I 

Plate 


Load  Cell 
(^102Xh25 


(Dimensions:  mm) 


Fig.  3  Dimensions  of  the  finger-crusher  test  rig 
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Table  1  Loading  conditions  for  the  component  loading  test 


Pulse  Duration, 
AT  (msec) 

Velocity  Change 
on  the  Loading  Table, 
AV  (m/s) 

Condition 

Code 

Velocity 

X04_30 

4 

4.5 

X04_45 

6.0 

X04 60 

3.0 

X08_30 

8 

4.5 

X08_45 

6.0 

X08 60 

3.0 

X20_30 

20 

4.5 

X20_45 

6.0 

X20 60 

2,3  Measurement 

The  arrangements  of  sensors  in  the  test  setup  are  shown  in  Fig.  4,  and  the  speeifieations  are  listed 
in  Table  2.  The  signals  from  sensors  were  reeorded  in  the  data  aequisition  system  (Speetral 
Dynamies,  SYSCHASVXI-5)  at  the  sampling  frequeney  of  250  kHz.  Moreover,  direet-eurrent 
offset  and  a  CFCIOOO  filter  were  applied  to  eaeh  signal  in  the  postproeessing  on  the  data- 
analysis  software  (MathWorks,  MATLAB).  The  test  was  also  reeorded  by  the  high-speed 
imaging  eamera  (Phantom,  Miro)  for  the  eonfirmation  of  the  speeimen’s  behavior  and  the  test 
rig’s  rigidity. 
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Potentiometer 
(Upper  Plate 
-Base  Plate) 


Loading  Table  / 

Load  Cell  /  Accelerometer  Accelerometer 

(Lower  Plate  (Loading  Table  A)  (Loading  Table  B) 

-  Base  Plate) 


Accelerometer 
(Lower  Plate) 


Fig.  4  Arrangements  of  sensors  in  the  test  setnp 


Table  2  Specifications  of  sensors  nsed  in  the  component  loading  test 


Sensor 

Position 

Type 

Range 

Accelerometer 

Loading  TableA 
Loading  TableB 

Endevco 

7270A-20KM6 

20000G 

Accelerometer 

Upper  Plate 

Lower  Plate 

Endevco 

7270A-60KM6 

60000G 

Load  cell 

Lower  Plate  -  Base  Plate 

Omegadyne 

LC8400-213 

lOOOOIbf 

Potentiometer 

Upper  Plate  -  Base  Plate 

Celesco 

MT2A-30E-33 

30in. 

The  input  acceleration  used  to  calculate  AV  on  the  loading  table  was  the  average  of  the  outputs 
from  2  accelerometers  (Loading  Tables  A  and  B).  The  load  applied  on  the  specimen  Lm  was 
calculated  from  the  output  from  the  load  cell  L  considering  the  inertial  effect  of  the  lower  plate, 
as  follows: 
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Lyn  =  L-miXai 


(1) 


where  mi  is  the  lower  plate’s  mass,  and  ai  is  the  aeeeleration  on  the  lower  plate. 


2,4  Finite  Element  Analysis 


There  are  3  FE-material  models  for  the  butyl  rubber  fabrieated  on  LS-DYNA  simulation 
software  by  different  developers:  Livermore  Software  Teehno logy  Corp.  (LSTC),  Humaneties, 
and  the  ARE.  In  this  study,  we  evaluated  the  aeeuraey  of  these  models  by  eomparing  their  FEA 
results  with  those  of  experiments  with  the  finger-erusher  test  rig.  Therefore,  the  FE  model  of  the 
eomponent  loading  test  was  built  on  LS-DYNA  (Fig.  5)  using  material  models  listed  in  Table  3 
for  the  test  rig’s  parts. 


Number  of  Elements:  58904  Mesh  Size  (Specimen):  2.5mm 


Fig.  5  FE  model  of  the  component  loading  test 


Table  3  FE-material  models  for  the  test  rig’s  parts 


Part 

Material 

Model  Type 
in  LS-DYNA 

Parameters 

Frame 

Upper  Plate 

Aluminum 

MAT_020_ 

RIGID 

RO=2.816e-6  (kg/mm^) 
E=70{GPa) 

Lower  Plate 

PR=0.30 

Load  Cell 

Stainless 

Steel 

MAT_020_ 

RIGID 

RO=8.000e-6  (kg/mm^) 
E=193(GPa) 

PR=0.30 
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Three  FE  material  models  listed  in  Table  4  were  used  for  the  speeimen  part  to  make 
eorresponding  FE  models.  The  ESTC  material  model  and  Humaneties  material  model  are 
derived  from  the  neek  eomponent  of  the  Hybrid  III  50*  ATD  FEA  models  developed  by  ESTC 
and  Humaneties.^’^  Both  models  are  derivative  of  a  viseoelastie  material  model,  and  the  ESTC 
model  is  the  most  simplified  deseription.^  In  addition  to  viseoelastieity,  the  Humaneties  model 
eontains  Ogden  hyperelastieity  for  the  equilibrium  part.'” 

Table  4  FE-material  models  for  the  butyl-rubber  specimen  part 


Developer 

Model  Type  in  LS-DYNA 

Model  Description 

LSTC 

MAT_006_ 

VISCOELASTIC 

Linear  Viscoelastic  Model 

Humaneties 

MAT_077O_ 

OGDEN_RUBBER 

Ogden  Hyperelastic 

Rubber  Model 
with  Viscoelasticity 
(First-O  rder  Prony  F  unction) 

ARL 

MAT_183_ 

SlvlPLFED_RUBBER 

_WrrH_DAMAGE 

Incompressible  Rubber  Model 
with  Strain-Rate  Dependent 
Loading  /  Unloading  Curves 

The  ES-DYNA  material  eards  used  for  the  ESTC  and  Humaneties  models  are  listed  in  Tables  5 
and  6.  By  eontrast,  the  ARE  model  is  a  form-type  model  and  refers  to  loading-unloading  eurves 
for  its  response  eonsidering  strain-rate  dependeney.^  The  material  eard  of  the  ARE  model  is 
listed  in  Table  7. 


Table  5  LS-DYNA  material  card  for  the  LSTC  model 


Parameters 

RO=1.1000e-6  (kg/mm^) 
BULK=0.1128(GPa) 
G0=0.0046  (GPa) 
Gl=0.0010  (GPa) 
BETA=0. 11000 
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Table  6  LS-DYNA  material  card  for  the  Humanetics  model 


Parameters 

RO=1.700e-6  (kg/mm^) 
PR=0.49 

MU1=8.100e-004  (GPa) 
MU2=7.220e-005  (GPa) 
MU3=-9.710e-005  (GPa) 
ALPHA1=1.3 
ALPHA2=4.0 
ALPHA3=-2.0 
GI=3.200E-2  (GPa) 
BETAI=1.600  (GPa) 


Table  7  LS-DYNA  material  card  for  the  ARL  model 

Parameters 

RO=1 .600e-6  (kg/mm^) 

K=2.900  (GPa) 

MU=0.900  (GPa) 

G=0.100  (GPa) 

SIGF=1.000e-4  (GPa) 


The  ARL  model’s  loading-unloading  eurves  are  developed  through  the  dynamie  loading  test 
developed  at  Purdue  University  as  shown  in  Fig.  6. 
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Fig.  6  Strain-rate-dependent  loading-nnloading  cnrves 
for  the  ARL  model 


Each  FE  model  was  applied  the  same  boundary  conditions  as  shown  in  Fig.  7  according  to  the 
acceleration  histories  on  the  loading  table  in  experiments.  Thus,  the  velocity  histories  in  Fig.  8- 
10  were  generated  on  FE  models  as  the  loading  conditions. 


Fig.  7  Schematic  of  honndary  conditions  in  FEAs 
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Fig.  8  Velocity  histories  on  the  loading  table  in  FEAs 
for  4-nullisecond  (msec)  pnlse  loading 
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Fig.  9  Velocity  histories  on  the  loading  table  in  FEAs 
for  8-msec-pnlse  loading 
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Fig.  10  Velocity  histories  on  the  loading  table  in  FEAs 
for  20-msec-pnlse  loading 


3.  Results  and  Discussion 


3,1  Experimental  Results 

Typical  load  and  deformation  histories  of  the  component  loading  test  are  shown  in  Fig.  11.  There 
is  a  delay  in  the  deformation  response  most  likely  due  to  the  viscoelastic  behavior  of  the  butyl 
rubber.  The  experimental  results  are  summarized  in  Table  8. 

The  deviation  of  the  results  in  each  loading  condition  is  less  than  10%  except  for  the  deformation 
peak  in  20-msec-pulse  loading,  and  these  are  believed  to  be  small  enough  to  compare  the  results 
with  those  of  FEAs.  Considering  the  accuracy  of  the  potentiometer  (Celesco  MT2A-30E-33, 
0.15%  of  full  scale),  it  seems  that  the  small  deformations  in  20-msec-pulse  loading  were  not 
measured  precisely.  Therefore,  the  reliability  of  these  results  is  lower  than  that  of  others.  This 
should  be  considered  in  the  following  discussion. 
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Fig.  11  Typical  load  and  deformation  histories  of  the  component  loading  test, 
X08_30_03  (AT  =  8  msec,  AV  =  3.0  meters  per  second  [m/s]) 
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Table  8  Summary  of  experimental  results 


Shot 

Code 

Pulse 
Duration, 
AT  (msec) 

Loading 
Table 
Velocity 
Change, 
AV  (m/s) 

Peak 

Load 

(kN) 

Applied 
Impulse 
(kN  msec) 

Peak 

Deformation 

(mm) 

X04_30_01 

4.7 

3.0 

13.3 

31.1 

3.36 

X04  30  02 

4.5 

3.0 

13.9 

32.1 

3.25 

X04_30_03 

4.2 

3.0 

13.4 

31.7 

3.14 

Avg.  /SD 

4.5/6% 

3.0/0% 

13.5/2% 

31.6/2% 

3.2513% 

X04_45_01 

3.9 

4.4 

22.7 

47.3 

4.87 

X04_45_02 

3.7 

4.4 

23.6 

47.6 

4.92 

X04_45_03 

3.9 

4.4 

22.6 

46.3 

5.09 

Avg.  /  SD 

3.8/3% 

4.4/0% 

23.0/2% 

47.1  /1% 

4.96/2% 

X04_60_01 

3.3 

6.0 

35.3 

62.4 

5.73 

X04  60  02 

3.3 

5.9 

36.2 

62.7 

6.22 

X04_60_03 

3.2 

5.9 

35.9 

63.2 

6.28 

Avg.  /  SD 

3.3/2% 

5.9/1% 

35.8/1% 

62.8/1% 

6.08/5% 

X08_30_01 

8.3 

2.9 

7.4 

29.7 

1.73 

X08  30  02 

8.5 

2.8 

7.2 

30.6 

1.93 

X08_30_03 

9.0 

3.0 

7.7 

30.4 

2.12 

Avg.  /  SD 

8.6/4% 

2.9/3% 

7.5  /  3% 

30.2/2% 

1 .93  /  9% 

X08_45_01 

8.4 

4.4 

11.3 

45.0 

2.77 

X08_45_02 

8.6 

4.4 

11.1 

44.3 

2.85 

X08_45_03 

8.0 

4.4 

11.2 

44.7 

2.66 

Avg.  /  SD 

8.3/4% 

4.4/0% 

11.2/1% 

44.7/1% 

2.76/3% 

X08_60_01 

7.9 

5.9 

16.1 

59.4 

4.01 

X08  60  02 

8.0 

5.9 

15.5 

58.6 

3.97 

X08_60_03 

8.0 

6.0 

16.1 

59.2 

3.74 

Avg.  /  SD 

8.0/1% 

5.9/1% 

15.9/2% 

59.0/1% 

3.91  /4% 

X20_30_01 

21.9 

3.5 

2.6 

32.3 

0.79 

X20_30_02 

21.7 

3.3 

2.6 

32.3 

0.63 

X20_30_03 

21.5 

3.3 

2.5 

30.0 

0.64 

Avg.  /  SD 

21.7/1% 

3.4/3% 

2.6/2% 

31.5/4% 

0.69/13% 

X20_45_01 

20.7 

4.6 

4.1 

44.2 

0.62 

X20_45_02 

19.9 

4.7 

4.1 

45.8 

0.60 

X20_45_03 

21.2 

4.8 

4.2 

45.4 

1.26 

Avg.  /  SD 

20.6/3% 

4.7/2% 

4.1/1% 

45.1  12% 

0.83/45% 

X20_60_01 

19.9 

6.5 

6.4 

62.0 

1.77 

X20  60  02 

19.3 

6.6 

6.8 

61.6 

1.81 

X20_60_03 

19.0 

6.5 

6.2 

60.6 

2.03 

Avg.  /  SD 

19.4/2% 

6.5/1% 

6.5/5% 

61.4/1% 

1 .87  /  7% 
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3,2  Comparison  in  Stress  and  Strain  Histories 

The  load  and  deformation  histories  of  eaeh  experiment  were  transferred  to  the  engineering  stress 
and  engineering  strain  histories  aeeording  to  the  measured  dimensions  of  the  speeimen.  The 
stress  and  strain  histories  are  shown  in  Figs.  12-17  with  eorresponding  FEA  results  of  eaeh 
material  model.  The  FEA  results  of  the  Humaneties  model  and  ARE  model  agree  with  the 
experimental  results  in  most  of  the  loading  eonditions  ineluding  delays  in  strain  histories. 


Time  (msec) 

Fig.  12  Example  of  the  engineering  stress  history 

comparison  between  experiments  and  FEAs, 
X04_60  (AT  =  4  msec,  AV  =  6.0  m/s) 
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Fig.  13  Example  of  the  engineering  stress  history 

comparison  between  experiments  and  FEAs, 

X08_60  (AT  =  8  msec,  AV  =  6.0  m/s) 

However,  the  FEA  results  of  the  LSTC  model  signifieantly  overshoot  the  experimental  results  in 
several  loading  eonditions.  Moreover,  the  LSTC  model  shows  unrealistie  double-eurved 
response  in  20-msee-pulse  loading  as  shown  in  Fig.  14  and  17. 


Time  (msec) 

Fig.  14  Example  of  the  engineering  stress  history 

comparison  between  experiments  and  FEAs, 
X20_60  (AT  =  20  msec,  AV  =  6.0  m/s) 
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Fig.  15  Example  of  the  engineering  strain  history  comparison 
between  experiments  and  FEAs,  X04_60 
(AT  =  4  msec,  AV  =  6.0  m/s) 


Time  (msec) 


Fig.  16  Example  of  the  engineering  strain  history  comparison 
between  experiments  and  FEAs,  X08_60 
(AT  =  8  msec,  AV  =  6.0  m/s) 
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Time  (msec) 

Fig.  17  Example  of  the  engineering  strain  history  comparison 
between  experiments  and  FEAs,  X20_60 
(AT  =  20  msec,  AV  =  6.0  m/s) 

To  clarify  the  accuracy  of  each  material  model,  the  errors  of  the  TEA  results  versus  the 
experimental  results  were  compared  for  the  peak  stress,  impulse  per  unit  cross-section,  and  peak 
strain  in  Figs.  18-20.  Thus,  the  errors  probably  have  some  relationships  with  the  pulse  durations 
because  the  tendencies  are  similar  among  loading  conditions  with  the  same  pulse  duration. 
Therefore,  we  will  discuss  the  accuracy  of  material  models  mainly  from  this  point  of  view. 

The  errors  of  the  Humanetics  model  and  the  ARE  model  are  less  than  25%  except  for  the  strain 
errors  in  20-msec-pulse  loading,  although  the  ESTC  model  has  considerably  larger  errors.  The 
strain  errors  in  20-msec-pulse  loading  are  still  large  after  considering  the  low  reliability  of  the 
deformation  measurement;  the  deviation  in  strain  histories  in  the  unloading  phase  of  20-msec- 
pulse  loading  is  also  large  (Fig.  17).  Because  of  the  characteristic  of  viscoelastic  material  models, 
the  Humanetics  model  assumes  the  same  behavior  for  the  loading-unloading  phases;  the  ARE 
model  considers  only  one  curve  for  the  unloading  phase  regardless  of  strain  rate  (as  was  shown 
in  Fig.  6).  Consequently,  the  unloading  phases  of  each  model  contain  larger  errors  than  those  of 
the  loading  phases  (Fig.  20).  Together  with  the  long  pulse  duration,  this  probably  leads  to  the 
accumulation  of  strain  errors  in  20-msec-pulse  loading. 

For  these  reasons,  both  the  Humanetics  model  and  ARE  model  can  be  used  for  the  prediction  of 
stress  and  strain  histories,  and  also  for  that  of  impulse  transfer  in  4-  and  8-msec-pulse  loading. 
However,  they  require  cautious  use  in  longer-duration  pulse  loading. 
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Fig.  18  Peak  stress  errors  of  FEAs  vs.  experimental  test  data 


Fig.  19  Impulse  per  unit  cross-section  errors  of  FEAs  vs.  experimental  test  data 
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Fig.  20  Peak  strain  errors  of  FEAs  vs.  experimental  test  data 
3.3  Comparison  in  Stress-Strain  Curves 

The  stress  and  strain  histories  of  eaeh  experiment  were  eombined  into  stress-strain  eurves.  The 
eurves  are  shown  in  Figs.  21-23  with  eorresponding  FEA  results  of  eaeh  material  model.  The 
stress-strain  eurves  in  the  experimental  results  show  different  paths  for  the  loading-unloading 
phases  and  form  eneireled  areas.  The  FEA  results  of  the  Humaneties  model  and  ARE  model 
agree  with  the  experimental  results  in  the  formation  of  these  eneireled  areas.  By  eontrast,  the 
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FEA  results  of  the  LSTC  model  show  almost  identieal  paths  for  the  loading-unloading  phases  in 
4-  and  8-msee-pulse  loading,  as  shown  in  Figs.  21  and  22. 


Engineering  Strain 


Fig.  21  Example  of  the  engineering  stress-strain  cnrve 
comparison  between  experiments  and  FEAs, 
X04_60  (AT  =  4  msec,  AV  =  6.0  m/s) 


Engineering  Strain 


Fig.  22  Example  of  the  engineering  stress-strain  cnrve 
comparison  between  experiments  and  FEAs, 
X08_60  (AT  =  8  msec,  AV  =  6.0  m/s) 
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Engineering  Strain 

Fig.  23  Example  of  the  engineering  stress-strain  cnrve 
comparison  between  experiments  and  FEAs, 

X20_60  (AT  =  20  msec,  AV  =  6.0  m/s) 

There  are  2  eases  for  the  formation  of  eneireled  areas  in  FEAs: 

1)  The  material  model  assumes  different  material  behaviors  for  the  loading-unloading  phases. 
However,  the  LSTC  model  is  a  viseoelastie  material  model  and  assumes  the  same  behavior. 

2)  The  material  model  assumes  different  material  behaviors  for  the  deformation  speeds,  and 
the  speeds  are  different  in  the  loading-unloading  phases.  This  is  the  ease  for  the  LSTC 
model,  and  will  be  diseussed  in  this  seetion. 

The  shear  relaxation  modulus  G(t)  is  defined  in  the  LSTC  model  as  follows: 

G(t)  =  Goo  +  (Go  -  Goo)  (2) 

where  Goois  the  infinite  shear  modulus,  Gq  is  the  short-time  shear  modulus,  and  is  the  deeay 
eonstant.  (Goo,  Go  and  P  eorrespond  to  the  parameters  GO,  Gl,  and  BETA  in  Table  5, 
respeetively.) 

Therefore,  the  sensitivity  of  the  material  behavior  against  deformation  speed  is  governed  by  in 
an  exponential  manner.  Judging  from  the  formation  of  the  eneireled  area  in  Fig.  23,  the  LSTC 
model  was  probably  tailored  for  the  loading  with  pulse  longer  than  20  msee.  Consequently,  the 
sensitivity  in  4-  and  8-msee-pulse  loading  is  not  enough  to  differentiate  the  loading-unloading 
eurves  aeeording  to  their  differenees  in  the  deformation  speeds. 

Although  the  Humaneties  model  has  a  similar  deseription  for  the  relaxation  modulus  as  the 
LSTC  model,  the  loading-unloading  eurves  are  differentiated.  Moreover,  the  ARL  model  also 
differentiates  the  loading-unloading  eurves.  Therefore,  the  Humaneties  model  and  ARL  model 
ean  be  used  for  the  evaluation  of  these  loading-unloading  events. 
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To  clarify  the  accuracy  of  each  material  model,  several  metrics,  Eo,  Eioad,  and  Sss,  are 
implemented  as  shown  in  Fig.  24.  Each  metric  is  defined  below. 


Stress,  a 


Fig.  24  Schematic  of  the  stress-strain  curve  and  metrics 


The  initial  modulus  Eo  represents  the  instantaneous  material  behavior  and  is  defined  as  follows: 


Eo 


gp-o 

0.005-0  ’ 


(5) 


where  ao  is  the  stress  at  0.5%  strain. 


The  loading  modulus  Eioad  represents  the  overall  loading  behavior  and  is  defined  as  follows: 


■•Load 


Ep-0 


(6) 


where  aMAx  is  the  maximum  stress,  and  sp  is  the  strain  at  the  maximum  stress  point. 


Encircled  area  Sss  represents  the  absorbed  energy  during  the  loading-unloading  event  and  is 
defined  as  follows: 

Sss  =  .  (7) 

where  cr’  is  the  stress  in  the  loading  curve,  a”  is  the  stress  in  the  unloading  curve,  and  smax  is 
the  maximum  strain. 


The  errors  of  the  FEA  results  versus  the  experimental  results  were  calculated  in  these  3  metrics, 
as  shown  in  Figs.  25-27.  The  errors  of  Eo  show  that  all  3  material  models  are  not  suitable  for  the 
prediction  of  instantaneous  material  behavior  such  as  the  rising  edges  of  the  experimental  results 
(in  Figs.  21  and  22). 
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■  LSTC  ■Humanetics  bARL  BExperimentSD 


Fig.  25  Initial  modulus  errors  of  FEAs  versus  experimental  test  data 


Fig.  26  Loading  modulus  errors  of  FEAs  versus  experimental  test  data 


Fig.  27  Encircled  area  errors  of  FEAs  versus  experimental  data.  (LSTC  model’s  errors  are  omitted  due 
to  poor  differentiation  in  the  loading-unloading  curves.) 

As  shown  in  Fig.  26,  the  errors  of  Eioad  are  less  than  25%  in  most  of  the  FEA  results  of  the  ARE 
model  in  4-  and  8-msee-pulse  loading.  The  only  exeeption  is  the  loading  eondition  X04_60;  this 
has  the  highest  strain  rate  due  to  its  short-duration,  high-amplitude  pulse.  The  strain  rates  of  the 
experimental  results  in  this  loading  eondition  are  shown  in  Fig.  28.  The  maximum  strain  rate  is 
approximately  150/s,  and  this  is  lower  than  the  upper  limit  of  the  strain  rate  defined  in  the 
loading  eurves  of  the  ARE  model  as  shown  in  Fig.  6  (160/s).  However,  the  strain  rate  in  Fig.  28 
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is  the  averaged  value  over  the  whole  speeimen.  Therefore,  the  partial  strain  rate  probably 
exeeeds  the  upper  limit,  and  the  ARL  model  tends  to  underestimate  as  Eo  in  this  loading 
eondition.  Moreover,  this  also  leads  to  the  underestimation  of  peak  stress  in  the  loading  eondition 
X04_60  as  shown  in  Fig.  18.  The  FEA  results  of  the  Humaneties  model  and  LSTC  model  also 
show  small  errors  of  Eioad  in  several  loading  eonditions,  but  the  errors  are  generally  larger  than 
those  of  the  ARL  model. 


Time  (msec) 

Fig.  28  Strain-rate  histories  of  the  specimen  in  the  test 
condition,  X04_60  (AT  =  4  msec,  AV  =  6.0  m/s) 

As  to  Sss,  the  Humaneties  model  and  ARL  model  show  small  errors  for  4-  and  8-msee-pulse 
loading.  However,  the  deviation  is  larger  than  that  of  Eo]  this  should  be  noted  in  the  evaluation 
of  energy  absorption  on  the  FEA.  In  20-msee-pulse  loading,  both  models  show  large  errors  due 
to  the  overestimated  strain  in  the  unloading  phase,  whieh  was  diseussed  previously. 

For  these  reasons,  both  the  Humaneties  model  and  ARL  model  ean  be  used  for  the  predietion  of 
stress-strain  eurves  in  4-  and  8-msee-pulse  loading,  whereas  the  instantaneous  material 
behaviors  are  not  eompletely  simulated.  Both  models  also  ean  be  used  for  the  predietion  of 
energy  absorption  with  larger  errors.  Moreover,  the  ARL  model  is  superior  to  the  Humaneties 
model,  espeeially  in  the  predietion  of  loading  behavior. 
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4.  Conclusions 


In  this  study,  we  conducted  component  loading  tests  of  butyl  rubber  and  compared  the  results 
with  those  of  FEAs  to  clarify  the  accuracy  of  the  FE-material  models.  The  Humanetics  model 
and  ARE  model  showed  good  accuracy  in  4-  and  8-msec-pulse  loading  for  the  prediction  of 
stress  history,  strain  history,  stress-strain  curve,  impulse  transfer,  and  energy  absorption,  while 
the  ESTC  model  showed  significant  errors.  Moreover,  the  ARE  model  is  superior  to  the 
Humanetics  model  especially  with  its  higher  accuracy  in  the  loading  phase.  Therefore,  we 
conclude  that  the  ARE  model  is  to  be  used  in  the  forthcoming  FEA  of  the  ATD’s  lower  leg  and 
the  butyl-rubber  floor  mat.  However,  the  current  ARE  model  has  two  problems  to  be  corrected: 

1)  The  strain-rate  dependency  in  the  unloading  phase  is  not  completely  installed  in  the 
models.  This  leads  to  the  excessive  recovery  delay,  especially  in  long-duration  pulse 
loading. 

2)  The  strain-rate  dependency  shown  in  Fig.  6  does  not  cover  the  partial,  high  strain-rate 
deformation.  This  leads  to  the  excessively  soft  behavior  in  the  loading  phase  and 
underestimation  of  the  maximum  stress. 

Problem  2,  in  particular,  affects  the  4 -msec-pulse  loading  that  we  frequently  use  in  the  evaluation 
of  underbody  blast-loading  events.  Therefore,  the  ARE  model’s  improvement  is  strongly 
recommended. 
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6.  List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

US  Army  Research  Eaboratory 

ATD 

anthropomorphic  test  device 

FE 

Finite  Element 

FEA 

Finite  Element  Analysis 

lED 

improvised  explosive  device 

kN 

kilonewton 

ESTC 

Eivermore  Software  Technology  Corporation 

mm 

millimeter 

msec 

millisecond 

m/s 

meters  per  second 
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1  DEFENSE  TECH  INFO  CTR 

(PDF)  ATTN  DTIC  OCA 

2  US  ARMY  RSRCH  LABORATORY 

(PDF)  ATTN  MAE  HRA  MAIL  &  RECORDS  MGMT 

ATTN  RDRL  CIO  EL  TECHL  LIB 

1  US  ARMY  RSRCH  LAB 

(PDF)  ATTN  RDRL  WMP  F  R  KARGUS 
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